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We report a density-functional study of the heteronuclear Coir dimer adsorbed on 
benzene or graphene. In either case Coir prefers an upright position above the center 
of a carbon hexagon with the Co atom next to it. The Ir atom stays away from the 
carbon ring and thus preserves its free-atom-like properties. This results in a very 
large magnetic anisotropy of more than 0.2 eV per dimer. So high a value should 
suffice for long-term data storage at the temperature of liquid nitrogen. 
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Graphene has recently attracted much attention in both applied and fundamental sci- 
ence. Owing to its novel properties, graphene opens new perspectives for post-silicon 
electronics^, spintronics^, chemical sensor-, hydrogen storage^, touch-screen panel^, and fur- 
ther applications^. Recently, our density-functional-theory (DFT) investigations revealed 
that both benzene and graphene are suitable substrates to adsorb some transition-metal 
dimers (TMD) in a magnetic state of potential interest for high-density magnetic data 
storage^. The hexagonal environment provided by the carbon ring preserves the electronic 
structure characteristics of the free Co dimer, which results in a magnetic anisotropy energy 
(MAE) of the order of 100 meV per dimer-. A systematic DFT study of a whole series of 
dimer-benzene complexes including the 3c? and 4c? TMD Fe2, C02, Ni 2 , RU2, Rh 2 , and Pd 2 
shows that promisingly large MAE, stable geometry, and stable magnetic ground state are 
also found in Ru 2 Bz (Bz = CgHg).- These findings may open a way to enhance the presently 
available area density of magnetic recording by up to three orders of magnitude. Consider- 
ing that a minimum MAE of 40 kT is required for long-term data stability^, the reported^ 
MAE could be sufficient to operate a magnetic storage device at about 30 K. The aim of the 
present work is to seek a possibility to further raise the MAE of TMD complexes in order 
to extend the expected operation range to the temperature of liquid nitrogen. 

Spin-orbit interaction in a magnetic state is the primary source of MAE^. Recent DFT 
calculations found that in the Co group, along with the increasing strength of the spin-orbit 
coupling parameter, the free dimer MAE increases from Co 2 to Ir 2 .— _ — Thus, depositing 
5c?-TMD like Ir 2 , on Bz or on graphene (Gr) might be expected to result in a larger MAE of 
the dimer-carbon system than in the related 3c?-TMD case. Previous studies however tell us 
that chemical binding to surfaces frequently reduces the magnetic moment of 4c? or 5c? metals 
due to their relatively small intra-atomic exchange (Stoner) integrals^. For example, the 
calculated ground states of Rh 2 or Pd 2 adsorbed on Bz are non-magnetic- Heterodimers 
of 3c? and 5c? elements would combine a large Stoner integral (3c?) with strong spin-orbit 
coupling (5c?). 

Electronic structure analysis suggests that only a perpendicular (upright) arrangement 
of the dimer on Bz or Gr (local symmetry Cq v ) is helpful to preserve the large MAE of the 
dimer-. In such a configuration, Fig. [T] (a), the elemental identity of the atom next to the 
carbon ring determines the binding characteristics, while the atom atop contributes mostly 
to the magnetic properties of the system. Therefore, choosing light atoms as base atoms and 
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placing heavy atoms atop may be a way to utilize the stronger spin-orbit coupling in the 
latter without sacrificing the high magnetic moment. Here, we use the heteronuclear dimer 
Coir to test this idea. The ground-state structures, binding energies, magnetic moments, 
and magnetic anisotropies of Coir adsorbed on benzene and graphene were investigated. 
For comparison and deeper understanding, these properties were also studied in Ir 2 Bz, in 
C02BZ, and in free Coir, Ir2, and C02 dimers. 

Our DFT calculations were performed with a highly accurate all-electron full-potential 
local-orbital scheme (FPLO)^, release 9.00-34^. The molecular mode of FPLO with free 
boundary conditions and Fermi temperature broadening of 100 K was used for free dimer and 
dimer-Bz systems. To simulate the interaction between the dimer and Gr, a 4x4 supercell 
(space group PI) of dimension 9.84x9. 84x16 A 3 was used, large enough to avoid interactions 
among dimers in neighboring cells. The related k-mesh used for linear-tetrahedron-method 
integrations contained 6x6x1 points in the Brillouin zone. All data reported here were 
obtained using the generalized gradient approximation (GGA) with the PBE96 exchange- 
correlation functional^. The valence basis set comprised Co (3s, 3p, 3d, 4s, 4p, Ad, 5s), Ir 
(4/, 5s, 5p, 5d, 6s, 6p, Qd, 7s), C (Is, 2s, 2p, 3s, 3p, 3d), and H (Is, 2s, 2p) states. Geometry 
optimization was carried out in a scalar relativistic mode with a force convergence threshold 
of 10~ 2 eV/A . 

To find the lowest-energy geometries and spin states of the CoIrBz complex, 18 kinds of 
initial structures were optimized without any symmetry constraints for total spin S = 0, 1, 
2, and 3. (See Fig. Al(i)-(vi) and (ix)-(xiv) of Ref . The structures (vii) and (viii) were 
ignored since there the distance between metal atoms is too short to be reasonable. For 
structures (i), (vi), (ix), (xi), (xii), and (xiii), 6 new variants appear by interchanging the 
positions of Co and Ir atoms.) Initial antiferromagnetic and ferrimagnetic states were also 
considered for S = and 1. 

For CoIrGr, 9 initial configurations were considered, denoted by || cc , || mm , || tt , -Ucoir, 
_L c i r Co, J-mCoir, J-mirCo, -Ucoir and -UirCo- Here, _L or II mean that the dimer axis is oriented 
perpendicular or parallel to the Gr plane, respectively; c, m, and t denote that Co and Ir 
are placed above the center of a carbon ring, above the midpoint of a C-C bond, or on top 
of a C atom, respectively. In the _L configurations, the atom mentioned first is next to the 
Gr plane. All carbon positions were fixed in the plane with a C-C bond length of 1.42 A. 

The MAE of CoIrBz or CoIrGr in the ground state was computed as the energy difference, 
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£[ioo] _ £[ooi] ^ between states with magnetization direction perpendicular [0 1] and parallel 
[1 0] to the carbon plane. The in-plane anisotropy was neglected, [1 0] was defined to 
point from the center to the corner of the carbon hexagon. In the case of free dimers, [0 1] 
and [1 0] stand for the directions parallel and perpendicular to the dimer axis. The MAE 
was evaluated in the fully relativistic mode in which spin-orbit coupling is included in all 
orders. As a matter of experience^, the MAE values obtained within standard GGA and 
by including the orbital polarization correction^ , give respectively a lower and an upper 
estimate of the expected MAE. 

Benzene is a good model system to describe the main characteristics of dimer adsorption 
on graphene^i. In Fig. [TJwe present the ground state and the four lowest-energy isomers of 
the CoIrBz complex, obtained by the described optimization. In the ground state, the Coir 
dimer is bound perpendicularly to benzene with Co next to the carbon ring (Cq v symmetry, 
Fig. Ufa)). The distance between Co and Ir in the free dimer is 2.13 A. This metal- metal 
bond length is increased by only 3% in the cases of upright adsorption (Fig. (Ha), (d)), 
while it is enlarged by about 10% in the parallel adsorption modes (Fig. Hfb), (c), (e)). The 
dimer adsorption energy, E a d, in the ground state is 1.51 eV, which is comparable to the E a< i 
value of C02BZ, 1.47 eV (Tab. I). We also notice that the calculated E a d of h^Bz, 1.74 eV, 
is larger than that of C02BZ or CoIrBz. For an upright adsorption with Ir next to Bz (not 
shown in Fig. [T|), we found a ferrimagnetic state with E ac i = 0.94 eV and S = 1 at almost 
the same energy as a ferromagnetic state with S—2. 

A detailed understanding of adhesion and isomer energies is complex and beyond the 
scope of this paper. It would require consideration of the transition-metal level positions in 
different spin states, their different orbital overlap and their element-specific Stoner integrals. 
To give just one example, we note that the Co contribution to the bonding 7r-states of CoIrBz 
situated at about -8 eV (Fig. [2]) is quite different in the two spin channels. This originates 
from the spin-splitting of the CoIr-7r states, right-hand panel of Fig. [21 and the resulting 
unequal energy differences from the Bz-HOMO level. 

The ground-state spin of CoIrBz, S = 2, equals that of the free Colr-dimer. This makes 
a welcome difference to Ir 2 , whose moment is completely quenched when it is adsorbed by 
benzene. Using the heteronuclear dimer Coir instead of Ir 2 protects the magnetism of the 
heavy atom since it stays atop and preserves its free-atom-like properties to a large extent. 
In the ground state of the CoIrBz complex, the magnetic moments of Co and Ir are about 
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1.8 fiB and 2.3 //#, respectively. Compared with the free Coir dimer, the magnetic moment 
of Co decreases while that of Ir increases (Tab. I). Fig. [2] shows the energy level scheme and 
the orbital composition of each level for the ground states of CoIrBz and of Coir. Like in 
C02BZ,- the adsorption is realized mainly by forming chemical bonds between carbon and 
the atom next to the benzene plane. For example, the 3d xy , 3d x 2_ y 2 orbitals of the Co atom 
and the LUMO of benzene form the 5 states at -5.9 eV/-1.7 eV (majority spin channel) and 
-4.8 eV/-1.5 eV (minority spin channel). The ir states around -8 eV in both spin channels 
mainly consist of the 3d xz , 3d yz orbitals of Co and the HOMO of benzene. At the Fermi 
level, a two-fold degenerate singly occupied 5* state is to 90% composed of the 5d xy , 5d x 2_ y 2 
states of Ir. It is the splitting of this state by the spin-orbit interaction that determines the 
MAE of the whole system. 

The value of MAE and the site-resolved spin and orbital moments fis an d for both 
magnetization orientations are listed in Tab. II. Data for free Coir and for CoIrGr (discussed 
below) are included for comparison. In all systems, the spin moments are nearly the same 
for [0 1] and for [1 0] orientation, but the anisotropy of the orbital moments is obvious. 
The Ir atom shows large orbital moments in the case of [0 1] orientation and relatively 
small values in the [1 0] case. The Co atoms behave differently in the free and in the bound 
dimers. In free Coir, the orbital moment of Co in the [0 1] orientation is around 0.9 
much larger than in the [1 0] orientation. When bound to benzene, the orbital moment 
of Co is quenched. This comes about through a re-distribution of Co-weight between the 5 
and 5* states in the minority spin channel, see Fig. [2J Interaction with Bz leads to a severe 
reduction of Co-weight in 5* in favor of the 5(5) and tt*(5) states. As a result, the singly 
occupied 5* molecular orbital, carrying a large ground-state orbital moment of nearly 2 
is situated almost exclusively on the Ir atom. 

The calculated lower and upper estimates of the MAE in CoIrBz are about 250 meV and 
290 meV, respectively. Both of them are larger than either of the values for free Coir. The 
reason is that the MAE is proportional to the orbital moment anisotropy and grows with the 
strength of spin-orbit coupling^ and that bonding to benzene results in a significant increase 
of the former on the Ir atom (Tab. II). 

Most notably, the lower estimate of the MAE of CoIrBz is five times higher than the 
respective estimate for C02BZ.- Because the spin-orbit coupling in Ir is much stronger than 
in Co, the MAE caused by the Ir orbital moment anisotropy in CoIrBz is correspondingly 
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much larger than the purely Co-related MAE in C02BZ. We also notice that C02BZ has 
a slightly higher upper MAE estimate than CoIrBz, which is due tc£ the larger Racah 
parameter B of Co (145 meV) than of Ir (94 meV). 

Finally, we investigated the properties of a Coir dimer deposited on graphene. Structural 
optimizations of all 9 initial configurations confirm that the Coir dimer is adsorbed by Gr 
in an upright position with local symmetry Cq v , and that the Co atom prefers to occupy the 
position next to the Gr plane. The calculated dimer adsorption energy for the ground state 
amounts to 0.63 eV, and the relative energy of two separate adatoms to a bound dimer, 
[-E-CoGr + -EirGr — E Gr ] — -^CoirGr, is 3.34 eV. This clearly indicates that the metal atoms prefer 
to aggregate as dimers rather than as separate adatoms. The Co-Ir bond length and the 
Co-C bond length are 2.18 A and 2.25 A, respectively, similar to the CoIrBz case. A band- 
structure analysis (not shown) for the ground state geometry and spin of CoIrGr indicates 
that the main contribution to the states near the Fermi level originates from the 5d x 2_ y 2 and 
5d xy orbitals of Ir. Like in CoIrBz, the Ir atom in CoIrGr shows free-atom-like properties 
and its orbital moment anisotropy gives rise to a very large MAE estimated to be between 
200 meV and 330 meV. 

In summary, the Coir hetero-dimer is expected to bind with benzene /graphene in per- 
pendicular orientation above the center of a carbon hexagon with the Ir atom farther away 
from carbon than the Co atom. Due to its remote position, the Ir atom largely preserves 
its free-atom-like properties. We predict that this kind of binding does not result in a de- 
terioration of the magnetic moment of the Coir dimer and that the MAE of this structure 
amounts to 0.2 . . . 0.3 eV. These numbers exceed all known experimental values of MAE by 
more than one order of magnitude. Judging by the criterion MAE > 40 CoIrGr should 
be a promising material for ultrahigh-density magnetic data storage at liquid-nitrogen tem- 
perature. The demonstrated similarity of our results for CoIrBz and CoIrGr suggests a 
possibility to use other substrates, like boron-nitride or various aromatic carbon-hydrates, 
with a related chance to achieve chemical self-assembling. 
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FIG. 1. (Color online) The dimer adsorption energy (E ad ), metal-metal and metal-carbon bond 
lengths (in A), and spin magnetic moments for the ground-state geometry (a) and higher energy 
geometries or spin states (b)-(e) of the CoIrBz complex. Hexagons indicate benzene rings. The 
red and blue bullets indicate Ir and Co atoms, respectively. 
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FIG. 2. (Color online) Scalar-relativistic single-particle energy levels and their orbital compositions 
for the ground states of CoIrBz (left panel) and Coir (right panel) . The upward (downward) arrow 
indicates majority (minority) spin states. All energies refer to a common vacuum level. Dimer- 
dominated states are labelled in black and benzene-dominated states are labelled in red. The three 
types of chemical bonds between the Coir dimer and Bz are labelled blue in parentheses. The 
position of the Fermi level (Ef) is indicated by a horizontal dashed line. For better resolution, the 
<5|-level of Coir is raised by 0.160 eV, the cr*-level is raised by 0.005 eV, and the 7T|-level is lowered 
by 0.012 eV. 
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TABLE I. Dimer adsorption energy E a d, total spin 5, atom-resolved spin magnetic moments, 
A*S(tmi) an d MS(TM2)> metal- metal and metal-carbon bond lengths, dTMi— TM2 and ^tmi-Cj for the 
ground state of C02, C02BZ, L"2, Ir2Bz, Coir, CoIrBz, and CoIrGr. TM1 denotes the atom next 
to the carbon plane and TM2 the atom farther away. The data for C02BZ differ slightly from our 
previous results^, because in this work full structural optimization was performed, while in Ref. Q 



the Bz plane was fixed. 
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TABLE II. Spin and orbital moments us and HL (in n B ) and MAE (E^ 1 - E^'P , in meV) of the 
ground states of Coir, CoIrBz and CoIrGr. Values calculated with and without the OP correction 



are given 


in the columns entitled SO+OP and SO, 


respectively. 
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